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Abstract 

We perform a linearization of the Poisson-Boltzmann (PB) density functional for spherical Wigner- 
Seitz cells that yields Debye-Hiickel-like equations agreeing asymptotically with the PB results in the 
weak-coupling (high-temperature) limit. Both the canonical (fixed number of microions) as well as the 
semi-grand-canonical (in contact with an infinite salt reservoir) cases are considered and discussed in a 
unified linearized framework. In the canonical case, for sufficiently large colloidal charges the linearized 
theory predicts the occurrence of a thermodynamical instability with an associated phase separation of 
the homogeneous suspension into dilute (gas) and dense (liquid) phases. In the semi-grand-canonical 
case it is predicted that the isothermal compressibility and the osmotic-pressure difference between the 
colloidal suspension and the salt reservoir become negative in the low-temperature, high-surface charge 
or infinite-dilution (of polyions) limits. As already pointed out in the literature for the latter case, these 
features are in disagreement with the exact nonlinear PB solution inside a Wigner-Seitz cell and are 
thus artifacts of the linearization. By using explicitly gauge-invariant forms of the electrostatic potential 
we show that these artifacts, although thermodynamically consistent with quadratic expansions of the 
nonlinear functional and osmotic pressure, may be traced back to the non-fulfillment of the underlying 
assumptions of the linearization. 

1 Introduction 

The study oLclassical charged systems has received an increased renewed interest in view of many industrial 
appUcationsiEfB paint, petrochemicals, food, pharmaceuticals, cosmetics, diapers, sewage treatment etc. 
Many environmental-friendly new materials are hydrosoluble due to the presence of ionizable groups that 
dissociate upon contact with water. In fact, their hydrosolubility is a result of the combination of Coulomb 
repulsion between fixed charged monomers and the mixing entropy maximized by the mobility in solution 
of the oppositely-charged small counterions. Besides technological applications, charged macromolecules 
like lipid aggregates (bilayers, micelles and vesicles), proteins and polynucleotides. lincluding DNA and 
RNA) are also of fundamental importance in the biochemistry of living systemsE'cl Furthermore, due 
to the availability of faster computers, many new insightaJn soft-matter physics come from Monte-Carlo 
and molecular-dynamics simulations of charged systems.BB These may be partially viewed as controlled 
theoretical experiments and provide a complementary approach to analytical treatments. 

An ubiquitous case is that of mesoscopic charged colloidal particles (also called polyions or macroions) 
immersed in aqueous solution, which polarize the small mobile ions in their vicinity: microions of opposite 
sign (counterions) are attracted to them, while like-sign microions (colons) are repelled. The theoretical 
description of these systems requires the understanding of the role of the electrostatic interactions between 
charged objects mediated by the surrounding aqueous ionic solution. In view of the many-body problem 
and the long-range nature of the Coulomb interaction, a statistical-mechanical treatment of the system is 
nontrivial. Within the Primitive ModelQ (PM) the molecular nature of the solvent is ignored (neglect of van 
der Waals and hydration forces) and the suspension is treated as a two-component system, comprised of the 
highly-charged large polyions (and its neutralizing counterions) and oppositely-charged pairs (anions and 
cations) of ionized salt particles. These are immersed in a continuous medium of dielectric constant e and 
interact through the bare Coulomb potential with additional hard-sphere repulsion. In the PM it is implicitly 
assumed that the (hard) spheres have the same dielectric constant as the solvent, so there are no electrostatic 
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image effects. For symmetric (in size and charge) electrolytes the PM reduces to the Restricted Primitive 
Model (RPM) and a theoretical descrixjlimi for dilute solutions may be developed using thelraditional Debye- 
Hiickel (DH) theory for electrolviesJj'QEHl with some improvements taking nonlinearitiealJ into account pi 
using integral-equation methods £3 An extension of these theories for a colloidal suspension is nontrivia]|13llj 
in view of the huge asymmetry between poly- and microions. Compared to the symmetric case, nonlinearities 
are magnified and dominate in the strong asymmetric colloidal limit. j— . j— . 

A mean-field approach to the PM, althcoicij-uoL-thermodynamically consistent Ji^'ta is represented by 
the Poisson-Boltzmann (PB) approximation.l!ii[liJL!-Jli3 This treatment gives a reasonable description in the 
weak-coupling (high temperature or small charge densities) limit, when the microionic correlations that 
are neglected at the mean- field level are unimportarit. But even the mean- field PB approximation for a 
suspension of charged polyions is a formidable taskcj due to the large asymmetry in size, mass and charge 
between the polyions and the small mobile ions. To jCircumvent this problem, the cell model has been 
introduced — originally for the cylindrical geometrycil' E3 — which reduces the many-body problem to the 
study of a (fixed) single polyion inside a Wigner-Seitz (WS) cell, whose volume is related to the polyion 
density. In the WS-cell model the single polyion plays only the role of a boundary condition. Note that 
this represents a dramatic simplification to the original PM formulation, where polyions and microions are 
treated at the same level. Even though the applicability of the WS-cell model appears at firstjSight to be only 
justified for an ordered crystalline phase, it has also been used to describe the fluid phase jHj'Eil neglecting 
thus the polyion translational entropy, the polyion-polyion interactions and intercell (both microion-polyion 
and microion-microion) correlations — intracell microion-microion correlations are neglected because of the 
mean-field approximation, which computes only the intracell microion-polyion correlations. To simplify 
further, the geometry of the WS cell is usually taken as to match the boundary condition on the polyion 
charged surface. For example, for colloidal particles a spherical charged polyion is considered inside a 
concentric spherical WS cell. This reduces the problem to the solution of a second-order ordinary (rather 
than partial) differential equation. For the spheiical geometry this-Leqiikes the numerical solution of the 
nonlinear PB equation, contrary to the planaro'tj'EII and cylindricaO'LJ'tS cases, when an exact analytical 
solution is possible. In analogy to the DH approach to the RPM, it would be thus very helpful to formulate 
a linearized version of the PB approximation for WS-cell models. We should remark, however, that the 
linearized version (at the mean-field level) of the WS-cell model does not include any intercell (neither 
polyion-microion nor microion-microion) correlations and intracell microion-microion correlations, contrary 
to the traditional DH approach to the (symmetric) RPM, which automatically includes them because in 
the RPM the mean-field contribU|tipn — which in the PB WS-cell model comes from the intracell polyion- 
microion correlations — vanishes.113 Therefore a more appropriate interpretation of the linearized equations 
to be obtained here is that they correspond to an expansion about the weak-coupling or high-temperature 
limit of the mean-field equations. 

However, expressions obtained within a linearized framework should be interpreted with caution, since 
they may lead to artifacts when applied outside their range of validity. Besides the specific case to be 
discussed in this paper, another example which is clearly an artifact of the linearization corresponds to the 
attractive component to the efiective interaction between two confined colloids induced by charged walls, prfiLi 
dieted under linearized theorynJ but in violation to the exact (at mean-field level) nonlinear PB repulsion. ^ 
Although earlier numerical analysis of the nonlinear solution were in agreement with the linearized theory]^ 
these were soon ruled out under very general conditions.EJ The disagreement with the rigorous nonlinear 
results might be attributed to flaws in the numerical calculations. —Attempts to include ionic correlations 
lead indeed to attractive contributions to the effective interactioixES'L^ but they are doubly-screened and 
thus are not able to overcome the repulsive electrostatic DLV0a'E3 (Derjaguin-Lapdau-Verwey-Overbeek) 
component. Therefore, experimentaLevidences of confinement-induced attractiont^EU'Ej and the occurrence 
of metastable piperheated crystalstJ can not be explained at the PB mean-field level and still remain an 
open question.Ej 

Almost two decades after the first experimental evidences of attraction between like-charged spherical 
colloids mediated by monovalent counterions in bulk deionized aqueous suspensions, its existence is still under 
dispute. Under the mentioned conditions, electrostatic-stabilized colloidal crystals have been investigated 
by Ise et al.l^ revealing the presence of empty regions (voids) inside the crystal. These experimental 
observations were interpreted as a coexistence between a dense crystalline phase and a dilute gas phase. 
Similar voids were also found experimentally in the fluid phaseJ13 which, in analogy to the critical behaviour 
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of symmetric electrolytes, were interpreted as a coexistence between dilute (gas) and dense (liquid) fluid 
phases. Even fully equilibrated macroscopic gas-liquid phase separation has been rep0rted,c3 although these 
experimental observations have been attributed to the presence of ionic impuritieso 

From the theoretical point-of-view attractive interactions between like-charged spheres are observed only 
under special conditions. jEcir |CX|ample, they have been seen in Monte Carlo simulations in the presence 
of multivalervUiCOunterionsE^Ea'EJ or when the low-temperature ordering of the discrete charges is taken 
into account .e2I Under the conditions described in the previous paragraph those controversial experimental 
findings are either aitributed to the presence of long-rangej-ataactive electrostatic interactions between 
like-charged polyionsea or by state-independent volume teiaasElp obtained by approximations that involve 
some kind of linearization: random-phase approximation,E3E3'xi DH pair-distribution functions-augmented 
by a variational approach for the polyicpipolyion interactions linear-response £mproximation,Ej extended 
DH thepry for asymmetric electrolytes)^ mean-spherical apprcpdmation (MSA)Ea and symmetpc PB and 
MSA.E2I Even though it has been argued by Overbeek and othersE3 that the Sogami-Ise attractionO is due to 
inconsistencies in thfiir thermodynamic treatment, the question does not seem to be settled yet and discussion 
is still in progress .O This_attractive potential is in contrast to the generally accepted repulsive electrostatic 
component of the DLVOEScJ pair potential between like-charged polyions. However, the purely repulsive 
nature of the polyion-polyion effective pair potepiial does not preclude a priori the existence of a liquid- 
gas separation, as has been shown by Roij et alL3 The focus on the polyion-polyion effective interactions 
overlooks the important contribution to the free energy due to the polyion-microion interactions. Because 
most of the alternative analytical calculations to the Sogami-Ise attractive interaction potential requires 
some linearization procedure, the predicted liquid-gas coexistence should be analyzed with caution. In 
these predictions disagree with simulation results in the presence of (explicit) monovalent counterions,! 
when no instabilities have been detected. Moreover, there are indications that the observed van der Waals- 
like loops are artifacts due to the linearization, these being lixastically suppressed when nonlinearities are 
reintroduced in the theory by the use of renormalized charges.Ej Furthermore the linearization of the WS-cell 
semi-grand-canonical PB functional — which describes at the mean-field level the system in electrochemical 
equilibrium with an infinite salt reservoir — yields negative-compressibility, thermodynamically unstable 
regions which are absent in a full nonlinear treatment.oE3 Although-|many aspects of these artifacts for 
the semi-grand-canonical case were alreadypjjeported in the literature,Efl including a general analysis of the 
linearization scheme for various geometries,Ea we believe that there are still a few subtle points that need to 
be clarified. 

The purpose of this paper is first to develop a linearization scheme suitable to the canonical (fixed amount 
of microions) case, by adopting an explicitly gautte.-invariant approach. For the semi-grand-canonical case, 
it has been shown by Deserno and von Griinberg^ that the occurrence of unstable linearized equations of 
state depends on the way the linearization scheme is performed and on the osmotic-pressure definition. By 
extending our gauge-invariant approach to the semi-grand-canonical ensemble, we try additionally to shed 
some light on this question. We argue that thermodynamic stability and consistency are in fact independent 
concepts. The gauge-invariant forms of the equations of state allow to establish an explicit correspondence 
between their nonlinear and linearized versions. The linearized equations, although thermodynamically 
self-consistent with quadratic expansions of the nonlinear ones, lead to artifacts when their underlying as- 
sumptions are not satisfied. We will show, by using gauge-invariant forms for the electrostatic potential, that 
there is a unique linearization (about the state-independent zero-th order Donnan densities) that corresponds 
to the minimization of the associated linearized semi-grand-canonical functional, which is also asymptotically 
exact (at the mean- field level) in the weak-coupling (high-temperature) limit. Therefore the expansion of 
the nonlinear functional about the state-independent Donnan densities, originally proposed for the spherical 
geometry by von Griinberg et alxi — and generalized for other geometries with analogous high symmetry in 
Ref. feS] — is not only "optimal" , but it is asymptotically exact in the weak-coupling limit. In a companion 
papercZI explicit analytical comparison is performed for the planar case, where the exact nonlinear solution 
(at the mean-field level) can be obtained. 

The remainder of the paper is organized as follows. In Section || the salt-free model is introduced and the 
associated nonlinear equations are reviewed. In Section ^ the linearization of the appropriate functional is 
performed, considering three distinct physical situations: the salt-free (in the presence of neutralizing counte- 
rions only) system introduced in Section ^ with fixed amount of added monovalent salt (canonical ensemble) 
and in electrochemical equilibrium with an infinite monovalent salt reservoir (semi-grand-canonical ensem- 



3 



ble) . In Section U we discuss our results in comparison to previous work and present some concluding 
remarks in Section ^. Appendix ^ presents the boundary-density theorem (at the nonlinear mean-field level) 
for the salt-free simplest case. In Appendix ^ it is shown that the linearized Helmholtz free energy may be 
obtained by a Debye charging process of the linearized electrostatic energy. Appendix |^ presents the formal 
derivation of the linearized osmotic pressure for the salt-free case, showing that it corresponds to a quadratic 
expansion of the corresponding nonlinear osmotic pressure. In Appendices ^ and ^ it is shown that the 
linearized osmotic pressure in the presence of salt for the canonical and semi-grand-canonical ensembles, 
respectively, are given by the quadratic expansions of the corresponding nonlinear osmotic pressures. In Ap- 
pendix ^ the self-consistent linearized averaged densities for the semi-grand-canonical ensemble are obtained 
by the minimization of the appropriate functional. In Appendix ^ we compare the linearized semi-grand- 
canonical osmotic pressures obtained by different schemes of the Legendre transformation connecting the 
canonical and semi-grand-canonical ensembles. 



2 Definition of the model 

Although_th£ derivation of the PB equation from a free-energy density functional can be found else- 
where J23'lla'E3 to introduce the notation and to stress the advantages of an explicitly gauge-invariant ap- 
proach, it is essential to rederive it in the following. For simplicity, in this section only the salt-free case (in 
the presence of neutralizing counterions only) is presented. Generalization of the model including monova- 
lent salt is straightforward and briefly described in Appendices ^ and ^ for the canonical and semi-grand- 
canonical cases, respectively. 

The system to be considered is a hard charged sphere (spherical polyion) of radius a and total charge 
—Zq distributed uniformly on its surface inside a concentric spherical WS cell of radius R > a, where q > 
is the elementary charge and Z ^ 1 is the polyion valence. The radius R of the WS cell is related to the 
polyion density Up such that the total volume of the suspension is equally distributed between the polyions, 
i.e. Up = (47r_R'^/3) ^ , whose hard cores occupy a fraction = (a/ R)'^ of the total volume. To ensure the 
overall WS-cell charge neutrality, there are Z mobile (positive) point-like counterions, each carrying a charge 
+9, that are allowed to move in the spherical shell a < |r| < i?, whose volume reads 

V= J d^r = ^ {R' - a') . (1) 

a<\r\<R 

Henceforth, to simplify the notation, it will be implicit that all integrations are performed over the free 
volume V unoccupied by the polyion core — but also including the charged surface at |r| = a — unless 
otherwise stated. At the mean-field PB level the counterions are treated as an inhomogeneous ideal gas and 
are described by their (continuous) average local number density n(r) = (I^J2i=i ^'^ ~ ''*^))' '^'^^^^ ^'^ 
three-dimensional Dirac delta unction and the double brackets denote an ensemble (Boltzmann-weighted) 
average over all positions of the counterions. The total charge number density (counterions plus the negative 
surface charge on the sphere), 

p{r)=n{r)-^S'i\r\-al (2) 

is related to the reduced electrostatic potential ^|J{r) = /3q'i>{r) by the (exact) Poisson equation, 

VV(r-) = -47r^Bp(r), (3) 

where £b = Pq'^/t is the Bjerrum length and /3^^ = fc^T is the thermal energy at temperature T. It is 
implicitly assumed that the solvent dielectric constant e remains the same inside the sphere, so image-charge 
effects due to dielectric contrast are absent. The formal solution to the Poisson equation (||) may be written 
in terms of the Green function G3{r,r'), 

V>(r) ^^B J Gsir, r') p{r'), V^Gsir, r') = -iirS^ir - r'), (4) 
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which in turn aUows us to express the mean- field Helmholtz free-energy functional ^[n{r)] of a single WS 
cell as 

(3.^[n{r)]= J dVn(r){ln[n(r)C^] + y J d^r d^r' p{r) G3{r,r') p{r'), (5) 

where C is the thermal de Broglie wavelength of the counterions. It should be remarked that the mean- field 
Helmholtz free-energy functional can be derived from the Jiiider lying PM Hamiltonian in different ways: 
as the saddle point of the action in a field-theoretic approach, Eaimm a Gibbs-Bogoljubov inequality applied 
to a trial product state that decouples the original HamiltonianE£l or fmm a first-cumulant expansion of the 
density functional reformulation of the associated partition functionO The first term of (||) represents the 
configurational entropy of the inhomogeneous ideal gas of counterions, while the second term corresponds 
to the electrostatic energy, which may be rewritten as 



^ J d3rdVp(r)G3(r,r')p(r') = ^J d^r^ir)p{r) 

d^r [\/tp{r)f - / 7/;(r) Wlr) • dS. (6) 

ottIb Jqv 



(3U 



The surface contributions to the electrostatic energy — the last term of (H), performed over the boundary 
dV of the free volume V — vanish due to Gauss' law and the overall WS-cell charge neutrality, 

Jd^rp{r)^0, or Jd^rn{r) = Z. (7) 

The equilibrium counterion-density profile is obtained by minimizing the mean-field functional J^[n(r)] 
with respect to n(r) under the constraint of overall WS-cell charge neutrality For this purpose it is 
convenient to introduce a translationally invariant (independent of r) Lagrange multiplier and to define 
the extended Helmholtz free-energy functional including a Lagrange-multiplier term. 



n[n{r)] = ^-p,i dV p{r) ^ ^ - p^i 



d^rn{r) - Z 



(8) 



which is the analogue of the conjugated semi-grand-canonical functional /2[n(r)] in the case of non-fixed 
number of counterions (but fixed number of polyions). We should not confuse the virtual Lagrange multiplier 
Pc\ with the a priori fixed chemical potential of counterions p in the presence of an infinite counterion 
reservoir, when the natural choice for the thermodynamic ensemble is the semi-grand-canonical one, 

n = ^ - p J d^rn{r), Jd^rn{r)^Z. (9) 

The Helmholtz free-energy functional ^ supplemented by the Lagrange-multiplier term, J7 7^ J7, does not 
correspond to the semi-grand-canonical functional J7, because J? and O are obtained under the constraints of 
a fixed number of counterions Z, and fixed chemical potential of counterions p, respectively. Although this 
distinction may be seem rather academic, it is important to stress that overall charge neutrality (ensured 
by the Lagrange multiplier pd) and electrochemical equilibrium (imposed by a fixed chemical potential of 
counterions p) are independent constraints. This prevents misinterpretations when treating the Donnan 
equilibrium (the system in electrochemical equilibrium with an infinite salt reservoir), when both constraints 
must be simultaneously satisfied. i— , 
Functional minimization of (3f2 with respect to n{r) leads to the Euler-Lagrange or stationary condition,E3 
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= ^(r) + In [n{r)C^] - (3p,, = 0, (10) 

n{r)—n{r) 



which yields the equilibrium counterion profile, n{r) = ePt^'^^-^i'^) j (^'^ , In the above formulas the equilibrium 
electrostatic potential, ip{r), is obtained by inserting the equilibrium counterion profile n{r) into ipir), 
Eq. (^ . The Boltzmann- weighted equilibrium counterion profile n(r ) , together with the Poisson equation (^) , 
yields the nonlinear PB equation, 

V2v;(r) = -l^c^^^-'^M + ^d\\r\ a), (11) 
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where the equihbrium Lagrange multipHer, 



/?/2ei = ln(neC')-ln(e"'^M), (12) 

is obtained by imposing the charge neutrality constraint (|^) . We introduced the effective average density of 
counterions in the free volume V unoccupied by the polyion core, 

n, ^ (n(r)> = |, (13) 

the brackets denoting unweighted spatial averages over V, 

, , f d^rX(r) , , 

^'^^^^^ " /d3r ' 

in contrast to their nominal mean density in the suspension, = nc(l — 0), which does not take into 
account the polyion hard cores. Note that the Lagrange multiplier (^2|) may be decomposed into two terms, 
a chemical potential corresponding to an ideal gas of uniform density ric and an electrostatic contribution due 
to the counterion-cloud polarization, which may be written in terms of an average related to the equilibrium 
electrostatic potential "^{r). We note that the presence of the Lagrange multiplier /Ici ensures explicitly 
the gauge invariance of the equations, that is, physical observables, like the equilibrium counterion-density 
profile, 

^e--^^"^ _ n,exp[(^)-V;(r)] 
/d3re-^M (exp[(V;)-V5(r)])' ^ ^ 

and the electric field, -E(r) = — V\E'(r) — — (/^g)"-^ VVj(r), clearly do not depend on a particular choice of the 
zero of the electrostatic potential '(/'(r), since observables can always be written in terms of the gauge- invariant 
difference /3/2ei ~ "^(t)- I'^ particular, explicitly gauge-invariant forms of the density profiles will be useful to 
establish a connection between the nonlinear and the linearized osmotic pressures (cf. Appendices ^ to ^ 
and to derive the self-consistent linearized averaged densities for the semi-grand-canonical ensemble (cf. Ap- 
pendix ^) . We should remark that — because consistent theories should always be gauge invariant — the 
use of explicitly gauge-invariant forms constitutes just a technical convenience, it does not represent an es- 
sential requirement. However, we believe that it provides a systematic and more transparent way to perform 
the calculations. Henceforth gauge-invariant will be a short-writing for explicitly gauge-invariant, unless 
otherwise stated. 

The mean-field Helmholtz free energy (of a single WS cell) F is obtained by inserting the equilibrium 
density profile n{r), Eq. ([l5|), into the mean- field functional ^, Eq. (|), 

PF = p,^[n{r)] ^ / d^r4>{r)n{r) - iz^(a) + f d^rn{r) {in [n(r)C^] - l} 



i / d^ri:{r)n{r). (16) 



In Appendix |a| it is shown that the (nonlinear) osmotic pressure P (over pure solvent), defined as the negative 
total derivative of the Helmholtz free energy F, Eq. ([l6|), with respect to the WS-cell free volumeEj V, is 
simply given by 



dV 



n{R), (17) 



Z,T 



which is the well-kiiow WS-cell mean-field result that the salt-free osmotic pressure is related to the counterion 
boundary densityE3 This simple-functional form still remains valid at the PM (beyond mean-field) level for 
WS-cells of various geometriesjid although the mean-field prediction for the equilibrium boundary density 
n{R) will (in general) disagree with the corresponding rigorous PM result due to the neglect of intracell 
microion-microion correlations and finite ionic sizes. Henceforth, to simplify the notation, we will omit the 
bar to denote equilibrium properties. 



6 



3 Linearization scheme 



3.1 In the presence of neutralizing counterions only 



Let us define a linearized free-energy functional that will lead to DH-like equations of state for the salt-free 
model system introduced in the previous section. We start by expanding the PB Helmholtz free-energy 
functional (||) about the uniform counterion density (^sj), i.e., up to second order in the difference n{r) — nc, 



/?^DH = Z [In (neC') - 1] 



1 d^r 


'n{r) ^ 




Tie 



The linear contribution is absent from (|l8|), because it automatically vanishes, 



d> 



n{r) 



- 1 



= V = 0. 



(18) 



(19) 



As in the previous section, we introduce a Lagrange multiplier /id and define the extended Helmholtz free- 
energy functional i7DH — ^dh — Moi J d^r p{r). Functional minimization of 17dh with respect to n{r) leads 
to the linearized equilibrium counterion profile, 



i{r) = nc [1 + (3fic\ - tp{r)] . 



(20) 



It follows from Eqs. ( |l9| ) and (|2^) that the linearized equilibrium Lagrange multiplier f3iici is given by the 
average electrostatic potential inside the WS cell. 



/3Mc1 = mr)) 



Inserting the equilibrium counterion-density profile (^0|) into the Poisson equation 
equation. 



(21) 

leads to the DH-like 



Zl 



V'^Pir) = [^{r) - /?//ei - 1] + ^S'{\r\ ~ a) 



(22) 



where the Debye screening length k ^ is evaluated with the uniform counterion density ric, 



K = Kc= \/iTr£Bnc. (23) 

This is different from the usual linearized-PB treatment of the spherical WS cell,!! where the Debye screening 
length is defined in terms of the WS-cell boundary density n{r = R). It should be remarked that the DH- 
like Eq. ( p^ ) leads to a gauge-invariant linearized solution, i.e., independent of the choice of the zero of the 
potential, because it depends only on the difference tp{r) — f3fic\- This is not the case for the commonly used 
linearized-PB solution for the spherical WS celljlil when it is disputable how to construct a self-consistent 
associated linearized free energy, from which the linearized equations may be obtained by a functional 
minimization. 

The solution of the linearized DH-like equation ( p^ ) for an arbitrary WS-cell boundary potential ip{R) 
and under the appropriate boundary conditions, 



dipir) 



reads 



dr 



ipir) 

f3qlBE{r) 

Ai(m,v) 
A2(m,v) 
A±(w) 



and 



d'0(r) 



dr 



= 0, 



r=R 



2Zk£b ZIb Ai{KR,Kr) 



A2{kR, na) r A2 (Ki?,Ka) 
_^ d^(r) ^ Z£lA2inR,Kr) 

dr r^ A2(Ki?, Kfl)' 

A+(M)e^- A_(w)e^", 

(l±u)e=F", 



(24) 

(25) 

(26) 

(27) 
(28) 
(29) 
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where the Lagrange multipHer PfJci is given by 



A2{kR, ko) 



(30) 



Inserting the DH-Uke solution for the electrostatic potential (|2^) into the linearized Helmholtz free-energy 
functional leads, after some algebra, to the linearized Helmholtz free energy, 



/3i^DH = /?^DH[«(r)] = Z 



In (neC') - - 



Z'^Ib Ai{KR,Ka) 
2a A2 {k R, Ka) 



(31) 



In Appendix ^ i±_is shown that the linearized Helmholtz free energy (|31|) may also be obtained by a Debye 
charging procesgli3 of the linearized electrostatic energy, confirming thus its thermodynamic self-consistency. 

The linearized osmotic pressure (over pure solvent) of the colloidal suspension follows from the negative 
total derivative of the linearized Helmholtz free energy (without including the Lagrange multiplier term) 
with respect to the WS-cell free volume V, 



dm 



DH 



DH 



Z.T 

Ai(kR, Ka) 



zHi 

2a 



d 

dV 



Ai(k_R, na) 



A2{kR, na) 



[Ai(Ki?,Ka) - A2{kR, 



Z.T 



Ana 1 



Zntf 



AA^{kR, na) 



' 1 2 
-K a 



K^R^ I - 



(32) 



where one should take into account both the explicit R dependence as well as the volume dependence 
of the screening length kT^ when computing the total derivative. The first term of (|3^ ) represents the 
uniform ccumterion-density ideal-gas law, while the next term corresponds to the mean-field electrostatic 
correctionaI3 due to intracell polyion-microion correlations. In Appendix ^ it is shown that the linearized 
osmotic pressure (^ ) can be also obtained by a formal differentiation of the linearized Helmholtz free energy 
and that it also corresponds to a quadratic expansion of the nonlinear PB osmotic pressure (0). At the 
end of the next subsection — that considers the system in the presence of added salt — we shall find that 
for sufficiently large surface charges the linearized osmotic pressure (^) is no longer a monotonic function 
of the WS-cell free volume V, which would imply a thermodynamical instability and an associated phase 
separation of the system — in contrast to the full nonlinear theory, which does not predict any instability. 



3.2 In the presence of neutralizing counterions and added salt (canonical en- 
semble) 

Let us now add a symmetric monovalent (1:1) salt to the system. At the same level of mean- field approx- 
imation, all microions are treated as inhomogeneous ideal gases, described by their average local number 
densities n±(r). We will not distinguish between counterions and positive ions derived from the salt dissoci- 
ation. Therefore n+(r) accounts both for counterions and positive salt ions (cations), while n_(r) represents 
the negative colons (anions). In terms of these number densities, the total charge number density reads 

p{r) = n+(r) - n_(r) - ^6'' (|r| - a) . (33) 

The effective average uniform densities of positive and negative microions are given by 

n± - (n±(r)) = Q+ = Z + nsV, Q- = n,V, (34) 

where rig is the a priori known effective average salt concentration and Q± are the fixed total number of 
positive and negative microions inside a WS cell. Within the cell-model approximation the salt ions are 
evenly distributed between different cells and the average salt concentration rig is the same for each identical 
WS cell. For later convenience, it will be useful to introduce the dimensionless parameter 



which measures the contribution of the sah ions to the ionic strength in the suspension, 

/ = i(nc + 2n,) = i(l + 2s)nc. 

For example, for s = 1 there is one pair of salt ions for each positive counterion in solution. 
We expand the nonlinear PB Hclmholtz free-energy functional, 

mn±{r)] = ^Jd'r [V^(r)]2 + J2 J dVn,(r) {in [n,(r)Cf] - l} , 
about the average uniform densities (^), up to second order in the differences n±{r) — h±, 
/3=^DH ("^C') - 1] + ^ / [V^(r)]^ + lY^n,Jd'r ^^^^^^ 



i=± 



(36) 



(37) 



(38) 



where C± are the thermal de Broglie wavelengths of cations (including the positive counterions) and anions, 
respectively. We introduce a Lagrange multiplier /ici, which will ensure the overall WS-cell charge neutrality. 



J d^rp{r) = 0, 



d^r [n+(r) - n_(r)] = Z, 



(39) 



and define the extended Helmholtz free-energy functional, 17dh — =^dh ~ /^ci J d'^r p{r). Functional mini- 
mization of J?DH with respect to n±(r) leads to the equilibrium density profiles. 



n±{r) = n± [1 ± /3/Zci T '^Pir)] 



(40) 



Using (H) and (|^), we obtain that the Lagrange multiplier is given by the average electrostatic potential 
inside the WS cell, /3/ici — (jp{r)), as in the salt-free case, Eq. (|2l|). Inserting the equilibrium profiles ( ^0[ ) 
into the Poisson equation (||), leads to the DH-like equation. 



■0(r) - Pfici 



1 



1 + 2s 

where the inverse of the Debye screening length is now given by 



Z£f 



+ ^SH\r\-a), 



(41) 



K = ^/SttIbI = V47r£B(l + 2s)nc = KcVI + 2s. 



(42) 



Solving the associated linearized DH-like equation leads again to the electrostatic potential (p5|), with the 
Lagrange multiplier for an arbitrary cell-boundary electrostatic potential ip{R) given by 



/3Mei = mr)) = V(i?) + 



2ZKiB 



A2inR,Ka) l + 2s' 



(43) 



At this point we should remark that — in the infinite-dilution limit (R — > oo, nc ^ 0), but in the presence of 
excess salt {s oo, finite ng) — the linearized electrostatic potential ipir) reduces to the Yukawa screened 
electrostatic potential. 



Z^rc-'^^''-'^) 



hm [V'(r) - V(i?)] = . , n = y^ri 



(44) 



which leads to the repulsive electrostatic component of the traditional DLVO interaction potentialBEl 
between two identical spherical particles of radius a whose centers are a distance r apart. 



W{r) = Z^li 



1 + KG 



(45) 



The phase diagram and dynamical properties of a suspension of spherical particles interacting through 
the effective DLVO pairwise potential (|5|) were systematically investigated in Ref. using molecular 
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and lattice dynamics techniques. We should note that the polyion-polyion interaction potential within the 
traditional (symmetric) DH framework may be obtained from the exact (non-spherically symmetric) solution 
of the Helmholtz equation V^ip{r) — K^-0(r) for two spherical charged particles.O The large-separation 
(r ^ oo) asymptotics of this pairwise potential leads directly to the DLVO interaction potential (|4^)..-Qne 
should also keep in mind that the exact limiting laws (within the RPM) of the underlying DH theorytSl — 
associated with the Helmholtz equation — are only valid in the asymptotic limit of vanishing ionic strengths 
(k 0). Most likely this exactness does not apply for the asymmetric case of strongly charged colloids. 
Alternatively, the DLVO interaction potential (^5|) may also be obtained by the large-separation asymptotics 
of the microion-averaged polyion-polyion potential of mean-force, oittained using the MSA integral equation 
for the polyion- and the microion-microion correlations in the PM.E3 

Inserting the equilibrium density profiles (MQ) into the linearized Helmholtz free-energy functional 



■^m[n±ir)], Eq. 



leads to the linearized Helmholtz free energy. 



/3i^DH = (3^DH[n±{r)] = (1 + s)Z {in [(1 + s)n,Cl] - l} + sZ [in {sn.C) - l] 



Z 
"2 



Z£b Ai{KR,Ka) 



A2iKR,Ka) l + 2s 



(46) 



The two first terms of (^6|) correspond to the ideal-gas entropy of the uniform expansion densities n± , while 
the last term represents the linearized excess Helmholtz free energy due to the polarization of the microionic 
cloud around the polyion. In the infinite-dilution limit and in the presence of excess salt {R oo, Uc 0, 
but finite Ug), 



lim — 



Z£b Ai{kR, Ka) 



1 



A2(Ki?, Kfl) 14- 2s 



2a (1 



BTLs 



(47) 



it coincides with the polyion-countcrion interaction free energy (including the polyion self-energy) obtained 
in an extended Debye-Hlickel-Bjerrum approach for colloidal suspensions — cf. Eq. (2) of Ref. |jl^. The 
asymptotic electrostatic excess free energy (^), obtained by linearization of the PB WS-cell model functional, 
accidentallyES coincides with the electrostatic work done in charging up the surface of the polyion against the 
ionic atmosphere in the framework of the traditionaLDH theory — c£_Ref. pp.339 — which is obtained 
by a Giintelberg charging process. The Giintelberglll and the DebyeE2l charging processes differ by the fact 
that in the latter the whole system (including the ionic atmosphere) is simultaneously charged. 

The connection between the infinite-dilution (i? — > oo) limit of the linearized excess Helmholtz free 
energy (ph. 



lim 



Z 
"2 



Zis Ai(k_R, Ka) 



1 



A2 (Ki?,Ka) l + 2s 



Z 
"2 



Zir 



1 



a(l + Ka) 1 + 2s 



(48) 



and the state-independent volume terms obtained by Roij et alL3 was first reported by Warren, cf. Eqs. (7) 
and (11) of Ref. [|B|, followed by Denton, cf. Eq. (55) of Ref. Q. Subtracting out from Eq. (^) the polyion 
self-energy, Z^^B/(2a), yields the negative contributions of the state-independent volume terms obtained by 
Roij et al. — cf. Eq. (61) of Ref. p3|| . 



ZIb 



1 



2 [a(l-fKa) l + 2s 



2a 



1 



— ^(Znp) , 



(49) 



recalling that n-p = (47ri?'^/3) is the polyion density of the suspension. It has been claimedEJ that 
these volume-dependent (but state-independent) negative contributions to the Helmholtz free energy are 
responsible for driving a gas-liquid phase separation in dilute deionized aqueous colloidal suspensions. 

The linearized canonical osmotic pressure of the colloidal suspension follows from the negative total 
derivative of the linearized Helmholtz free energy (without including the Lagrange multiplier term) with 
respect to the WS-cell free volume V , but keeping fixed the total amount of salt. 



PP^^i^.s) 



dm 



DH 



= (1 



dV 
- 2s)n, 



Z,T,s 



(1 + 2s)nc - 
Zk£b 



zHi 

2a 



_d_ 

dV 



Ai(k_R, Ka) 



4(H-2s)A2(Ki?, Ka) 



A2{kR, Ka) 
Ai(kR, Ka) 



Z,T,s 



Ka 
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X [Ai{kR, no) - A2{kR, kg)] - Ana ( 1 + " 



(50) 



In Appendix ^ it is shown that the hnearized canonical osmotic pressure 
expansion of the nonhnear PB canonical osmotic pressure (D5). 

In the vanishing volume fraction of polyions (infinite-dilution) limit, (j) 
canonical osmotic pressure has the asymptotic behaviour 



corresponds to a quadratic 
(a/R)^ — > 0, the linearized 



(3Pi 



can 
DH 



'(l + 2g)0 



10(1 + 2s) 



4^ 
175^ 



which leads to the asymptotic linearized inverse isothermal compressibility, 



1 

DH 



d/3R 



can 
DH 



= ZTip (1 + 2s) 



29 



15(1 + 2s) 



/,i/3 



4^ 
105' 



2/3 



e{<t>) 



(51) 



(52) 



In contrast to the semi-grand-canonical case (to be treated in the next subsection) , the linearized isothermal 
compressibility in the canonical ensemble is stable in the infinite-dilution limit, lim^^o Xdh > 0. However, 
as shown in Figures ^ and ^, for finite densities (0 ^ 0) and sufficiently large values of 0, the linearized 
canonical osmotic pressure (3P^^ is no longer a convex function of the volume fraction 0, implying thus the 
onset of a thermodynamical instability. For salt-free suspensions (s = 0), the associated critical point — 
represented by the black circle in Figure — is located at 



0.008586189- •■ 



44.902477094 • ■ 



(53) 



which is determined by the criticality condition dP^/d(f> = d^P^/d(lP' = 0. We should stress that the 
coexistence regions between the dilute gas (G) and the dense liquid (L) phases — not shown in Figure || — 
must be determined under the constraints of constant chemical potential of polyions and of salt particles 
Us- For example, the salt-free bimodal line that defines the boundary of the coexistence region is given by 
the coupled system of equations. 



^5h( 



6g,s = 0) 



^dh( 



^'l,s = 0), 



^lp{(t)G,s = 0) = ^J.p{cj)h, s = 0), 



where fip{(f), s = 0) is the salt-free chemical potential of the polyions. 



^p(0,s = 0) = 



DH 



V 



/),s = 0) + 



s = 0) 



s=0 



(54) 



(55) 



and V = 47ra'^/(30) = l/up is the total WS-cell volume. Because the critical behaviour is a spurious result 
of the linearization, it is not worthwhile to construct the phase diagrams in detail and we restrict ourselves 
only to present the spinodal lines in Figure ^. This, however, will not change the location of the critical 
points, where the bimodal and the spinodal lines meet. Addition of salt stabilizes the suspension by reducing 
the unstable region to higher values of 0, as also illustrated by the s ^ spinodal lines in Figure |^. The 
corresponding bimodal lines are determined by generalizing Eqs. (|53) to the case of added salt (s ^ 0), 



can / 
DH I 



.,s), 



Aip('/'G,s) = Aip(0L,s), 



(56) 



where the chemical potential of the polyions in presence of added salt, /ip((/), s ^ 0), is given by Eq. (G6). 
With added salt, there is (in the coexistence region) a redistribution of microion pairs between the dilute 
and the dense phases, which still obeys the WS-cell charge- neutrality constraint (p9[). Hence, inside the 



coexistence region {(f>Q < (j) < 0l) the two fluid phases will not have the same value of s (corresponding to 
the homogeneous system). Their values in the gas (sq) and the liquid (sl) phases are obtained by imposing 
the total conservation of salt particles. 



(0l - 0g) 



(t>G < (t> < 4>l- 



(57) 



and the equality in both fluid phases of the chemical potential of salt particles /is(^, s) — as given by Eq. (G5), 

Ms(0g,sg) = Ms(0l,sl). (58) 
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An important and useful concept in charged colloidal systems is the charge renormalizatiorOEj of the 
polyion bare charge Z for finite volume fractions.ci'. For highly charged polyions, i.e. in the o = Z/(47ra^) — > 
cxD limit, it has been shown by Alexander et aZ.Ej that the renormalized effective charge Zotin_lhe salt-free 
system approaches a saturation value Zgat ~ clw{4')/^b, with assuming numerical valuesoEj around 9 to 
15 in the volume-fraction range 0.01 < (j) < 0.1. A self-consistent linearized osmotic pressure including charge 
renormalization effects would require the inclusion of additional terms due to the volume- fraction dependence 
of the effective charge Z^ff = Zcs(V), since the osmotic pressure is defined as the negative total derivative 
of the Helmholtz free energy with respect to the volume V. In other words, the linearized canonical osmotic 
pressure taking into account charge-renormalization effects is not simply given by replacing Z Zcff 
into Eq. (|50|). This point will be considered in a future workEl It is interesting to note, however, that the 
linearized critical threshold ^crit given by Eq. ( |53| ) is very close to the (laxgest) salt-free saturatic|tt.(Z oo) 
effective charge ^sat = 3.Zsat^B/a ~ 45 determined by Alexander et alc£ It has been speculatedll3 that this 
curious coincidence drives suspensions of highly charged colloids close to criticality. This might acccpiil_fpi; 
some of the experimental findings in dilute deionized aqueous suspensions of highly charged collaids,!!!!^!!^ 
which then would be explained by the presence of strong density fluctuations near the criticalityJij Another 
consequence that effective charges are below the saturation value, and therefore also below the linearized 
critical threshold, is that charge renormalization would stabilize the suspension against phase separation, 
because the unstable region predicted by linearized theory is unreachable (or at least^drastically reduced) 
when including renormalized effective charges. This fact was pointed out previouslyEI using a generalized 
Debye-Hiickel-Bjerrum approach for charged colloidal suspensions. In the present calculation, however, 
the critical behaviour is an artifact of the linearization, which is absent in the full nonlinear treatment. 
Therefore, the occurrence of a thermodynamical instability can only be explained beyond the WS-cell mean- 
field approximation, by including excluded-volume effects, intercell polyion-microion and microion-microion 
correlations that are neglected in the WS-cell mean-field PB picture. 



3.3 In contact with an infinite salt reservoir (semi-grand-canonical ensemble) 

Let us now consider the colloidal suspension in electrochemical equilibrium with an infinite salt reservoir of 
fixed bulk density Uh- The suspension is separated from the infinite reservoir by a semi-permeable membrane. 
The solvent and microions (counterions and salt ions) can pass through the membrane, but not the large 
polyions. This gives rise to an imbalance in the osmotic pressure across the semi-permeable membraae—jt 
equilibrium between the suspension and the salt reservoir is referred to as a Donnan equilibrium H^eSc 
Like in the previous subsections we will consider only the case of monovalent counterions and symmetric 
monovalent (1:1) salt. 

The effective average salt concentration in the colloidal suspension, Ug = (n_(r)), does not coincide 
with the reservoir bulk density rib and is not known a priori. A nontrivial question is its dependence with 
the physical parameters of the system, e.g., bulk salt concentration n^, polyion radius a, polyion valence 
Z and volume fraction ip = {a/R)^. At thei-WS-cell PB mean- field level of approximation this problem 
has already been considered in the literatureE3EZl and it is summarized in Appendix ^. Compared to the 
canonical case treated in the previous subsection — when the amount of microions is fixed and known a 
priori — there are two main differences. First, because the Donnan equilibrium is established under constant 
electrochemical potential of microions, the natural thermodynamical ensemble to perform the calculations 
is the semi-grand-canonical one, 

/2dh = ^DH ~^lJ-t d^r n,(r), /3^± = In (nbC±) , (59) 

i=± •' 

where we impose the equality of the microion electrochemical potentials inside the colloidal suspension, ^± , 
to the (mean-field) chemical potential of ideal gases of uniform density rib for both types of ions in the 
infinite salt reservoir, In (nbC±)- The second difference is that the effective average uniform densities of 
positive and negative ions, about which the linearization should be performed, vary in a nontrivial way as 
the WS-cell free volume V is changed. In Appendix ^ it is shown that the self-consistent linearized average 
densities for the Donnan problem are given by the state-independent zero-th order Donnan densities, 



^nl + (2nb)^ ± 

n± = ;^ , (60) 
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where ric is the effective averaged countcrion density, defined in (^^. These correspond to the uniform 
densities that the system would have in the infinite-temperature {is — 0) hmit under the constraint of 
overall WS-cell charge neutrality (|39|). We should remark that they do not correspond to the effective 



averages of the full nonlinear PB densities (E3), 



{n±(r)) — 2 — 2 ' 

because of the nonvanishing quadratic and higher-order {v > 2) contributions of the electrostatic potential 
differences, 

S,ir)^m-il;{r)r . (62) 
Here one can see an advantage of the gauge-invariant formulation of the problem: the linearized expan- 



sion densities (60) are simply obtained by taking the potential-independent, infinite-temperature {£b = 
0) limit of the PB nonlinear averaged densities (|6^), when the electrostatic potential differences vanish, 
lim^g_+o (5;y(r) 0. It is also clear that any other choice for the expansion densities will not lead to the 
exact potential-independent (infinite-temperature) limit of the nonlinear equations. As also derived in Ap- 
pendix ^, the Donnan average densities, Eqs. (|60|), do indeed correspond to the functional minimization 
the corresponding linearized semi-grand-canonical functional. We note that Deserno and von Griinberj 
define an optimal linearization point V'opt — related to the above gauge-invariant expansion densities by 
fij. = nbe^'^°p* — by using arguments based on the plausibility of this choice. They show that any other 
choice for the linearization point would lead to confiicting inequalities involving nonlinear and linearized 
averages. Finally, we should mention that even though the uniform expansion densities (|6^) are internally 
self-consistent (within the semi-grand-canonical ensemble) with the linearized truncation of the quadratic 
expansion of the nonlinear semi-grand-canonical functional, global self-consistency (between the canonical 
and the semi-grand-canonical ensembles) will require the use of the quadratic truncation of the nonlinear 
averages (|l]) as the self-consistent averaged densities, as discussed in detail in Appendix The inclusion 
of the quadratic contributions into the expansion averages, however, do not improve the agreement between 
the linearized and noiilinear equations, as can be shown by the explicit analytical comparison in the exactly 
solvable planar case.cZl 

Once again, we introduce a Lagrange multiplier /id that enforces the overall WS-cell charge-neutrality ( ^9[) , 
and define the extended semi-grand-canonical functional, SIuh — ■!?dh— Mci / d^rp(r). Functional minimiza- 
tion of i?DH with respect to n±{r), performed in Appendix]^, leads to the self-consistent linearized averaged 
densities (|60| ) and the linearized equilibrium profiles, 

n±(r) = n±[l±(^A(r))T^(r)]. (63) 

We should remark that, contrary to the salt-free and canonical cases, cf. Eq. (^l|), here the equilibrium 
Lagrange multiplier /3/Ltei is not related to the average electrostatic potential, /3/ioi 7^ (V^('"))- 

Inserting the linearized equilibrium profiles (^^ into the Poisson equation (^), leads to the DH-like 
equation, 

V^i^ir) = ^2 [^(^) „ ^^(^)^ _ ^] ^ -^S^lrl ~ a), (64) 

where the parameter 

"-++"-- -I- (2nb)2 

measures the relative importance of the counterions to the ionic strength in the colloidal suspension, 

I^Un++ n-) = i + (2nb)2 = ^ = — (66) 

Furthermore, the (effective) Debye screening length k^^ in the colloidal suspension, 

2 2 

= SttIbI ='^= , '^^ > kI, (67) 
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is always shorter than the Debye screening length = l/v^STrFgnb associated with the bulk salt concen- 
tration rib in the reservoir, showing that screening is enhanced in the colloidal suspension compared to the 
salt reservoir. 

Solving the associated linearized DH-like equation (|6j) leads again to the electrostatic potential (p5|), 
with the average electrostatic potential over the WS-cell volume for an arbitrary cell-boundary electrostatic 
potential V-'(^) given by 

(^(r)) = + , , - V (68) 



A2{k,R, ko) 

In close analogy to the linearized Helmholtz free energy for the canonical case, Eq. 



3) , the final expres- 



sion for the linearized semi-grand-canonical potential is 
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(69) 



The linearized semi- grand- canonical osmotic pressure follows from the negative total derivative of the 
linearized semi-grand-canonical potential (without including the Lagrange multiplier term associated with 
/ici) with respect to the WS-cell free volume V, but keeping fixed the microion chemical potentials fi± = 



DH 



dV 
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(70) 



where the prefactor of the right-hand side represents the ideal-gas Donnan osmotic pressure taking only the 
WS-cell charge neutrality (p9|) into account, 



V 



(2nb) 



(71) 



In analogy to the salt-free, Eq. (|3^), and the canonical, Eq. (^), cases, the first term of ( [70|) represents 
the ideal-gas law associated to the state-independent zero-th order Donnan densities, while the remaining 
terms correspond to the mean-field electrostatic corrections due to the microionic polarization around the 
polyion. We should note that the second term inside curly brackets depends only on 77 and is thus 
independent. This could suggest that the zero-th order Donnan ideal-gas law, lim^g^o /3^dh ~ "c/?7, would 
not be recovered in the weak-coupling limit. This i?B-indcpendcnt term, however, is indeed necessary to 
cancel the contributions that arise from the last term in the Ib —* ^ limit in order to give the correct 
potential-independent infinite-temperature limit. Moreover, it is shown in Appendix ^ that the linearized 
semi-grand-canonical osmotic pressure ( [70| ) corresponds to a quadratic expansion of the nonlinear semi- 
grand-canonical osmotic pressure (E5). 

The linearized osmotic-pressure difference between the colloidal suspension and the infinite salt reservoir 

('T-c/??)\/l 



obeys /3APdh = /3-Pdh " = /3-Pdh " Wh)\J^ - ^f^ with PP^^ given by Eq. (^. In the next section 
we show that the linearized semi-grand-canonical osmotic-pressure difference /3APdh is intrinsically thermo- 
dynamically unstable in thej_infinite-dilution limit and we compare with expressions previously obtained by 
Deserno and von Griinberg.tS 



4 Comparison with Deserno and von Griinberg results 

As already pointed out previously in the hterature0'@ in this section we will show that the linearized 
semi-grand-canonical osmotic pressure defined by Eq. (|7^) yields artifacts in the low-temperature, high- 
surface charge or infinite-dilution (of polyions) limits. In contradiction to the exact nonlinear result ( |E5| ), 
which yields an osmotic-pressure difference that is always positivejES j3AP = (3P — 2nb > 0, the linearized 
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version APdh becomes negative in the above mentioned limits. In an attempt to define the osmotic pressure 
in a linearized framework, Deserno and von GriinbergtJ introduced an alternative definition, Pi, that has the 
advantage of being exempt from any instabilities and is obtained via the partial derivative of the semi-grand- 
canonical potential with respect to the volume, keeping the optimal linearization point ■(/'opt (to be defined 
below) fixed. Their second osmotic pressure definition coincides with the linearized version (|70| ) obtained 
in the previous section, P2 = P^^, and is obtained via the total derivative of the semi-grand-canonical 
potential with respect to the volume. These two distinct osmotic-pressure definitions are given, respectively, 
by Eqs. (43) and (44) of Ref. pi for d = 3, 
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(72) 
(73) 



2cosh^Vopt I 60 

V = h/2{'^R) Ky^i^a) - Ky^ii^R) 4/2 («;a), (74) 
£ = h/2{>iR) Ky^i'^a) + Ky^inR) hM'^a), (75) 

where {I^,K^} are the modified Bessel functions^ of the first and the second kind, respectively, and the 
optimal linearization point "iAopt satisfies the relations 
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In accordance with Eqs. (23) and (26) of Ref. pq] , they can be recast in a simpler formal form in terms of 
the gauge-invariant electrostatic potential differences Siy{r), defined by Eq. (p2h, 
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(77) 
(78) 



from which one can see that they differ by a term that is quadratic in the electrostatic-potential difference. 
Looking at Eq. (C14) of Appendix one may trace back that the omitted contribution in Pi originates from 
the volume dependence of tius optimal linearization point i/'opt, in accordance to the interpretation given by 
Deserno and von GriinbergE3 for the two distinct pressure definitions. We should recall that the linearized 



semi-grand-canonical osmotic-pressure ([70|) coincides with the second pressure definition, Pqj^ = as 
shown in Appendix ^ by a quadratic expansion of the nonlinear osmotic pressure. It corresponds indeed to 
the negative total derivative of the linearized semi-grand-canonical potential Odh with respect to the WS-cell 
free volume V, which we thus believe to be the consistent and correct definition of the osmotic pressure. 
It is convenient to introduce the dimensionless linearized osmotic-pressure differences, 



2nb 2nb 



i = 1,2. 



(79) 



In the vanishing volume fraction of polyions (infinite-dilution) limit, ~ {a/Kf' — > 0, we may write the 
asymptotic linearized osmotic-pressure differences in terms oi 9 = 'iZlBja, and a = Kbfl, 



Hi - 

n2 = 



0^6|202/3exp (-2a0-i/3) 



56* 



2a^ 



5 (1-1- a)' 



- 1 



,6*2 02/3 (^_2a(i 



-1/3 



which lead to the asymptotic linearized inverse isothermal compressibilities. 
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(83) 

While neglecting the contribution of the last quadratic term in the linearized osmotic pressure P\ always leads 
to positive isothermal compressibilities, lim^^o Xi > 0; its inclusion in Pi always yields negative isothermal 
compressibilities in the infinite-dilution limit for nonvanishing is-, liiii0_>o X2 < 0. This means that the 
pressure definition P2 predicts that the infinite-dilution phase is unstable, in contrast to the canonical case, 
as shown in Figure ^. Therefore the thermodynamically consistent linearized osmotic pressure P2 intrinsically 
fails at infinite dilution, leading to negative isothermal compressibilities in this limit. This fact was first 
noticed for the spherical case in Ref. and generalized to WS cells in d-dimensions in Ref. For 
finite densities, in contrast to the canonical case (Figure |2|), the low-cf) dilute (gas) phase may only be stable 
for sufficiently small bulk salt concentrations, as shown in Figure ^. For sufficiently large the finite- 
temperature critical point — located in the vicinity of the salt-free critical point (the black circle in Figure ||) 
— disappears. Note the nonmonotonic behaviour of the spinodal lines associated to P2 for KbO = 10~^ and 
Kbfl = 10^ (inset), which leads to oscillating osmotic pressures. We should note, however, that the WS-cell 
model ceases to be meaningful for such high volume fractions in the latter case of strong screening. Finally we 
should remark that beyond the linearized PB WS-cell model approximation the thermodynamical instability 
at the infinite-dilution limit will be removed by taking into account the translational entropy of the polyions, 
which yields an osmotic-pressure contribution that is linear in (j> and, therefore, overcomes the negative cubic 
leading term in the asymptotic linearized osmotic- pressure difference (|8l|). However, because our analysis 
restricts to the linearization of the PB WS-cell model, the effect of this stabilizing entropic contribution — 
which may drastically alter the spinodal lines, specially in the low-volume fraction region — is not considered 
here. Note that the PB WS-cell model, in its full nonlinear version, is fully stableEj even without invoking 
this stabilizing contribution. 

Let us stress again that thermodynamic consistency and stability are independent concepts. This can be 
illustrated by inspecting the two linearized osmotic-pressure definitions proposed in Ref. | |65t . The linearized 
osmotic pressure 112, although not fully stable, is self-consistent with quadratic expansions of the nonlinear 
osmotic pressure. The unstable region of 112 just reflects the breakdown of the linearization scheme to the 
nonlinear PB equation. On the other hand. Hi, that displays positive isothermal compressibilities in the 
infinite-dilution limit and does not present any instabilities, is not thermodynamically consistent, its stability 
being purely fortuitous. This can be further fortified by comparing the exact analytical expressions of the 
nonlinear osmotic pressure for the planar geometrycZI with the two corresponding linearized versions. In the 
subsequent paper we show that both linearized expressions. Hi and 112, approach asymptotically the exact 
nonlinear result in the appropriate (weak-coupling, £b — + 0) limit. Their convergence, however, are very 
different and it is the self-consistent definition 112 that gives a better approximation to the full nonlinear 
equation. Although analytical proofs can only be obtained for the planar case, we believe the same argument 
applies for any geometry. 



5 Concluding remarks 

We performed a linearization scheme consistent with quadratic expansions of the appropriate nonlinear 
thermodynamic functional. By using gauge-invariant forms of the electrostatic potential, we have shown 
that the linearized osmotic pressures correspond to quadratic expansions of the corresponding nonlinear 
versions for the three cases investigated: in the presence of neutralizing counterions only (salt-free case), 
in the presence of fixed amount of added salt (canonical case) and in electrochemical equilibrium with an 
infinite salt-reservoir (Donnan eawiiibrium, semi-grand-canonical case) . 

Contrary to previous works,ESEa we adopted a gauge-invariant formulation with the inclusion of a La- 
grange multiplier term to account for the charge-neutrality constraint. In the case of the Donnan equilibrium, 
it is shown that the minimization of the associated linearized semi-grand-canonical functional leads indeed 
to the state-independent zero-th order Donnan densities as the self-consistent expansion averages for the 
linearization. Therefore the optimality of the optimal expansion point ipopt introduced by Deserno and von 
GriinbergES can be understood as corresponding to a self-consistent minimization of the linearized functional. 



16 



It is shown that the self-consistent Hiifiarized osmotic pressure in the semi-grand-canonical ensemble — 
as already pointed out in the literatureOu — leads to artifacts in the infinite-dilution, high-surface charge 
and strong-coupling limits, where it predicts negative osmotic-pressure differences between the colloidal 
suspension and the salt reservoir and negative isothermal compressibilities. Attempts to define a fully stable 
linearized equation of state-lDased on the partial derivative of the linearized semi-grand-canonical potential 
with respect to the volumeE3 can not be justified in our approach based on the minimization of the linearized 
functional, its stability being a fortuitous result. This can be seen most clearly in the analytically tractable 
case of two infinite charged planes in electrochemical equilibrium with an infinite salt reservoir, which we 
present in a companion paper E3 

To avoid confusion we should stress at this point the exactness of the PB nonlinear solution at the 
mean- field level, its range of validity and limitations. In this work we discussed the linearization procedure 
in the framework of the nonlinear PB and the WS-cell model. The linearization constitutes here an ap- 
proximation to the nonlinear treatment, whose exact results (at mean- field level) may be then compared 
to the linearized ones, allowing a control over the approximations and the onset of possible artifacts in- 
troduced by the linearization. Of course, we are not able to predict correct results for real systems when 
the (starting) nonlinear theory itself breaks down. In this case, an eventual linearized result may acci- 
dentally lead to a correct prediction of, say, phase separation, due to the simultaneous application of two 
inadequate approximations, namely, the mean-field PB equation and its subsequent linearization. The fact 
that PB nonlinear theory for WS-cell models always leads to stable suspensions does not invalidate phase 
separation in real systems, which may be due to finite-size effects, intra- and intercell microion-microion 
and intercell polyion-microion correlations that are neglected in the WS-cell mean-field PB picture. In the 
Donnan-equilibrium case, one should also take the microion-microion correlations in the infinite salt reservoir 
into account. Hnwpvpr flip ph']."^ spparatinn in dciouized aqueoiiiSjiiiispensions of charged colloids predicted 
by linearized JihpirieaEj'OOEj'Ej'QBEj — which are claimedEj'Ej to be driven by the state-independent 
volume terma£3Ej — seems to be related to mathematical artifacts of the linearizatinn itself and does not 
correspond thus to a real physical effect, as already pointed out in previous workaHj'u via the linearization 
of the PB functional for WS-cell models and revisited it here. 
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A Nonlinear osmotic-pressure boundary density theorem 

In this Appendix it will be shown that the salt-free nonlinear osmotic pressure is related to the counterion 
density at the WS-cell boundary — see also Section 3 of Ref. [|l9| for the fixed-gauge derivation. For the 
Donnan problem, the fixed-gauge derivation is presented in Section II. C of Ref. ||65| . 

The osmotic pressure (over the pure solvent) P is defined as the negative total derivative of the Helmholtz 
free energy F — the functional ^ evaluated at the equilibrium profile n{r) = h{r) — with respect to the 
WS-cell free volume V, 

dF _ d^[n(r)] 

The calculation of the osmotic pressure starting from the explicit form of the optimized Helmholtz free 
energy F, Eq. (^6|), turns out to be nontrivial, because both the equilibrium counterion profile fi(r) as 
well as the Lagrange multiplier /3/ici depend implicitly on the WS-cell free volume V. However, because 
/9/iei was introduced as a Lagrange multiplier to enforce the charge neutrality, it is much easier to consider 
the total derivative with respect to the volume of the extended Helmholtz free-energy functional i?[n(r)] = 
^ — /id / d"^r p(r) — but now considering the Lagrange multiplier /id independent of the WS-cell free volume 
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V — evaluated at the optimized profile n{r) = n{r), 



dV 



sv 



df}[n{r)] 



i(r) 



dV 



i(r) 



^3^, 5n[n{r)] 



(5n(r') 



l{r') 



dn(r') 



(A2) 



where the derivatives of O are taken with constant /ioi- Because J?[rt(r)] is stationary with respect to the 
optimized profile n{r), 



6n(r') 



0, 



(A3) 



i(r') 



only the partial-derivative term contributes to the osmotic pressure. Furthermore, because the only explicit 
dependence of /?[7i(r)] comes from the integration limit,ll3 we are lead to 



(3P 



dV 



df3n[nir)] 



dV 



dp^[n{r)] 



n{r) 



dv 



n(r) 



n(r) 



= -^j{R)n{R) - n{R) {in [n{R)C^] - 1} + Pfld n{R) = n{R), 



(A4) 



which is the mean- field resultB that relates the salt-free nonlinear osmotic pressure to the counterion density 
at the WS-cell boundary when intracell microion-microion correlations are neglected. 



B Linearized free energy from a Debye charging process 

In this Appiendix it is shown that the linearized Helmholtz free energy can also be obtained by a Debye charg- 
ing procesalj of the linearized electrostatic energy. This is obtained by inserting the DH-like solution ( p5[ 
into (|), 



i3Uin,eB) 



d^r [Vxp{r)Y 



2aA^{KR,Ka) 



Ai{kR, Ka)A2{KR, Ko) — K^a / drA^{KR,Kr) 

J a 



(Bl) 



The linearized excess Helmholtz free aiergy over the ideal-gas entropy of the reference uniform state is then 
obtained by a Debye charging processt2l of the linearized electrostatic energy (Bl), 



/3i^DH - Z [In (neC') - 1] - 2 



A 



Z^£b 
2a 



Ai{kR, ko) 



3a 



A2{KR,Ka) K^{R^-a^) 



(B2) 



which regains Eq. (|3l|), obtained by adding the linearized electrostatic energy (Bl) to the quadratic expansion 
of the ideal-gas entropy. 



C Formal derivation of the Unearized osmotic pressure 

In this Appendix it will be shown that the salt-free linearized osmotic pressure is given by a quadratic 
expansion of the nonlinear counterion density at the WS-cell boundary. 



d/3FDH dp^uHHr)] 



dV 



dV 



ne |l + (^) - m) + \ [{i>) - ^(i?)]' - \ ( K'^) - i^X 



(CI) 



which may be obtained by truncating the expansion of the salt-free nonlinear PB osmotic pressure up to the 
quadratic terms. 



(C2) 
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where the bar denotes equihbrium properties. One route to check that the salt-free hnearized osmotic 
pressure ( |3^) may be cast in the form ( |Cl[ ) is by using the hnearized DH-hke solution to the electrostatic 
potential ( P5[ ) and computing explicitly the spatial averages 



Ai{KR,Kr)\ _ An 
Al{KR,Kr)\ _ An ''^ 

V 



/ drr Ai(/s:_R, Acr) 

J a 
rR 
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An A2{kR, nr) 
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A2{k.R, na) 



2ZIf 



Ai{kR, no) 



X [Ai(kR, ko) - A2(Ki?, Ko)] - Ana (l - k^R^) - Ak^R^ 



(C3) 



(C4) 



Alternatively, Eq. ( plD may be checked by taking a formal functional derivative with respect to the 
WS-cell free volume V — analogously as performed in Appendix ^ for the nonlinear treatment — of the 
salt- free linearized Helmholtz free energy functional (|l^), which is recast in the form 



/3=^DH[n(r)] = Z [in(neC') - l] + j d?r f[n{r)], 

n{r) 



(C5) 
(C6) 



where p{r) and ipif) defined in Section ^jby Eqs. (^) and respectively. The equilibrium counterion 
profile n{r) that minimizes the linearized version of the extended Helmholtz free-energy functional i?DH = 
=^DH - Mci / d^rn(r), with Pfl^i = {-ip), is given by 



n{r) = nc [l + (V') - i/^ir)] . 
As explained in Appendix |^ for the nonlinear PB treatment, it is much easier to evaluate 
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dV 



because the second term vanishes due to the stationary condition. The first term reads 

dBilY)H\n(r)] ^ d f , , d f duc f ,q df\n(r)] 
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Using the explicit form ( p7| ) of the linearized equilibrium counterion profile n(r), we obtain 
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(C9) 
(CIO) 
(Cll) 

(C12) 
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(C14) 



which yields the salt-free linearized osmotic pressure Eq. (|Cl| 
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D Linearized osmotic pressure: canonical ensemble 



In this Appendix it will be shown that the linearized canonical osmotic pressure is given by a quadratic 
expansion of the nonlinear canonical osmotic pressure. 
The nonlinear PB Helmholtz free-energy functional, 

/3^[n±(r)] = ^ J d'r [V^{r)f + J2 J dVn,(r) {in [n,(r)Cf] - l} , (Dl) 

is minimized with respect to n±{r) under the constraint of overall WS-cell charge neutrality, Eq. (^). 
We introduce a Lagrange multiplier /id and define the extended Helmholtz free-energy functional, /? = 
^ — /ici /d'^rp(r). Functional minimization of /? with respect to n±{r) yields the Boltzmann factors, 
n±{r) = e±^^'"T'^(''VC±- The Lagrange multiplier, 

/3/iei = In {n+CD - In (e-'^'^'')) = In (n_C-) + In (c^^'')) , (D2) 

is obtained by using the effective average definitions (|3^), which automatically satisfy the overall WS-cell 
charge neutrality, Eq. (^9|). Replacing into the Boltzmann factors, we obtain the nonlinear equilibrium 
density profiles, 

^ n± e'FV'('-) ^ n± exp [± (ip) T V'(t-)] , . 

"^^""^ (eT^W) (exp [± (V-) T ^(r-)]) ■ ^ ' 



The simple ideal-gas relation (|17|) between the nonlinear osmotic pressure and the WS-cell boundary density 
can be generalized in the presence of added salt. In this case, it can be shown that the nonlinear osmotic 
pressure is given by the total microionic density, n = + n_ , evaluated at the WS-cell boundary. 



dBF 



= n{R) ^ + + p4 

z,T,s (e-'/'f'')) (e'/'('')) 



dV 

Let us now compare the linearized canonical osmotic pressure (|50|) with a quadratic expansion about the 



uniform reference state of its nonlinear counterpart, Eq. (D4) 



pF = [1 + s)nc 



■ (exp m - V;(r)]) (exp [- (V-) + V(r)]) 

= (1 + 2s)n, |l + Y^-^ili?) + \52{R) - \ {52{r)) + G MR), {Ss{r))]^ , (D5) 

where the electrostatic potential differences S^{r) are given by (|6^ ). Because of the redefinition of the 
screening length k~-^ in terms of k~^, the spatial averages ( |C3| ) and ( |C^ needed to evaluate and 
will be multiplied by a factor (1 -I- 2s)~^. Using the explicit DH-like solution to the electrostatic 
potential (p5|), it is indeed possible t o sh ow that the linearized canonical osmotic pressure ( |50| ) corresponds 



to the truncation of the expansion (D5) up to the quadratic terms, i.e. (3P = PP^m + ^['^3(^)1 ('^3('''))]7 
with /JP^ff — »t.dh(^) + (1 + 2s) (nc/2) [62{R) — (<52(^))], where the first term corresponds to the sum of the 
linearized canonical densities at the WS-cell boundary. We see that the linearized canonical osmotic pressure 
is not simply given by the linearized boundary density riDH(P)7 because of the presence of the quadratic 
terms. 



E Linearized osmotic pressure: semi-grand-canonical ensemble 

In this Appendix it will be shown that the linearized semi-grand-canonical osmotic pressure is given by a 
quadratic expansion of the nonlinear semi-grand-canonical osmotic pressure. 

The nonlinear PB semi-grand-canonical functional associated to fixed microion chemical potentials (3ii± = 
In (nbC±), 
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^ J d'r [Wi,{r)f+J2 J d'rn,{r) jln 



ni{r) 



rib 



(El) 



is minimized with respect to n±{r) under the constraint of overall WS-cell charge neutrality, Eq. i^^. 
We introduce a Lagrange multiplier /Zd and define the extended semi-grand-canonical functional, i? = 
^ ^ Mci / d^r p(r). Functional minimization of i? with respect to n±{r) yields the Boltzmann factors, 
n±{r) = ribe^^^"'^'^'-''^. The Lagrange multiplier. 



+ (2nb)2 (e'^W) (e-V-W) ± 
2nb(eT'/'W) 



is obtained by imposing the overall WS-cell charge neutrality, Eq. 
factors, we obtain the nonlinear equilibrium density profiles. 
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Replacing into the Boltzmann 
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where the electrostatic potential differences (5^(r) are given by ( 

simple ideal-gas relation between the osmotic pressure and the total microionic density, n 
WS-cell boundary is still valid, leading to 



At the nonlinear mean-field level the 
n+ + n-, at the 
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Let us again compare the linearized semi-grand-canonical osmotic pressure (|70| ) with a quadratic expansion 
about the zero-th order Donnan densities n±, Eqs . (|60|), of its nonlinear counterpart, Eq. (EE). Now, because 
of the redefinition of k, the spatial averages (|C3| ) and (C4) needed to evaluate (ip) and {ip'^) are multiplied 
by a factor r/. Using again the explicit DH solution (^), it is possible to show that the linearized semi-grand- 
canonical osmotic pressure (|7^) corresponds to the truncation of the expansion ( |E5| ) up to the quadratic 
terms, i.e. PP = /SP^'g + MR), {S^ir))], with /JP^'g - n^R) + [nj{2rj)] [<52(P) - {62^)] , where 
the first term represents the sum of the linearized semi-grand-canonical densities at the WS-cell boundary. 
We see again that the linearized semi-grand-canonical osmotic pressure does not correspond to the linearized 
boundary density nj^n^R), because of the presence of the quadratic terms. 



F Linearized averaged densities for the semi-grand-canonical en- 
semble 



In this Appendix it will be shown that the self-consistent linearized averaged densities for the semi-grand- 
canonical ensemble are given by the state-independent zero-th order Donnan densities, which are obtained 
by a self-consistent minimization of the linearized semi-grand-canonical functional under the WS-cell charge- 
neutrality constraint. 

To obtain the self-consistent averaged densities up to second order we need to expand the nonlinear 
semi-grand-canonical functional, Eq. (El), up to third order about the (a priori unknown) effective average 
densities {n±{r)) = (l/V) J d'^rn±{r), 



f30[n±{r)] = ^ J d'r [V^(r)]Vy^i 



In- 



+ ^[(n,)ln (n,)] Jd^rS,{r) 



J dV<5,^(r)-i^(n,) J dV <5,3(r.) ^ [| d'rSf{r) 



(Fl) 
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where we introduced the relative deviations about the averaged densities, 



o±(rj = -, - 1. 



(F2) 



Keeping only quadratic terms in Eq. ( p^ ) is similar in spirit to the quadratic expansion presented in Eq. (13) 
of Rcf. |65| . However, we want to stress that, in general, the unknown average expansion densities {n±) 
depend itself on the ionic profiles. Functional minimization of Eq. with respect to the profiles n±{r) 
must take this fact into account, in addition to the WS-cell charge-neutrality constraint, Eq. (|3^). Eventually, 
for a linearized theory, the self-consistent expansion densities n± turn out to be indeed independent on the 



ionic profiles, cf. Eq. (60), but this can only be derived a posteriori. 



We introduce a Lagrange multiplier /Xei to ensure the overall WS-cell charge neutrality (|39|), and define 
the extended semi-grand-canonical functional, J? = 17 — /id J d^r p{r). Functional minimization of i7 with 
respect to n±{r) leads to the Euler-Lagrange or stationary conditions, SO /6n±{r) = 0, which may be cast 
in the form 



n±{r) = (n±) { 1 T [V'(r) - /Smci] - In 



1 



4(r-) + ^[4W,(4W)] 



(F3) 



where we have neglected the cubic averaged contribution, ^ [(<5±(t))] , because the neglected quartic term 
of (Fl) will also contributed to it. To obtain the averages (n±) self-consistently, we integrate n±{r) over the 
volume to obtain the consistency relations. 



(n±(r)) = exp \ T [(V^(r)) - M + ^ (SUr)) + 6 [(4 W)] 



(F4) 



where the Lagrange multiplier /ici is found by imposing the overall WS-cell charge neutrality (p9|). We 
should stress that in addition to the chemical potential of microions /i that defines the semi-grand-canonical 
ensemble, we introduced a Lagrange multiplier fid, whose role is twofold: besides the overall WS-cell charge 
neutrality (^), it also ensures the gauge invariance of the electrostatic potential ip{r). We should not confuse 
the chemical potential of microions /i, which is fixed by the bulk salt concentration rib of the reservoir, with 
the linearized Lagrange multiplier /id, which is associated with the Donnan effect and ensures the charge 
neutrality of the WS cell ( ^ ) . 

Before we derive the averages consistent with a linearization of the PB equation, let us also obtain the 
self-consistent averages corresponding to a quadratic approximation of the nonlinear equation. Noting that 
^±{1") — ±Si{r) + ff[S2{r), {S2{r))] and neglecting cubic terms in Eqs. (F2) and (F4) leads to the quadratic 
self-consistent averages, 

v/n? + (2nb)2e(*2W)±nc 



(n±(r)) 

and the quadratic equilibrium density profiles. 



n±(r) = (n±(r)) 



l±5i(r) + i52(r)-i(<52(r)) 



(F5) 



(F6) 



where the electrostatic potential i^-th order difference 6^{r) is defined by (|62[). These correspond indeed to 
the quadratic expansions of the nonlinear PB average densities and equilibrium density profiles, respectively, 
and are correct up to S2ir) (quadratic) terms. Here we may see another advantage of the gauge-invariant 
formulation: it provides us a systematic way to consider self-consistent higher-order approximations of the 
nonlinear equations, while the fixed-gauge analysis of Deserno and von GriinbergLJ does not allow them to 
extend their calculations to include higher-order terms. 

However, in order to be consistent with a linearization of the PB equation one needs to neglect also 
the quadratic terms in the approximate Euler-Lagrange conditions (F3), although global self-consistency 
of the Legendre transformation will require to include them — cf. Appendix |^. It is clear that ignoring 
these terms will yield the state-independent zero-th order Donnan densities as the self-consistent linearized 
averaged densities, 

V'n^ + (2nb)2 ± ric 



(n±(r)) 



(F7) 
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and the linearized equilibrium density profiles, 



n±{r) = {n±{r)} [1 ± Si{r) 



(F8) 



In the main text we used the notation h± = (n±)]^, Eq. (60), to refer to the linearized self-consistent 
averaged densities for the semi-grand-canonical ensemble, where the subscript '1' emphasizes the fact that 
the average densities were obtained under linearization. Deserno and von Griinberg^ justify this choice for 
the expansion densities — written in terms of an optimal linearization point "tpopt defined by n± = nhe^^°^* 
— by arguing that any other choice for the linearization point would lead to conflicting inequalities involving 
nonlinear and linearized averages. In a gauge-invariant formulation, however, the justification is indeed based 
on the self-consistent minimization of the linearized semi-grand-canonical functional ^dh[^±(^)], which is 
obtained by truncating the expansion of the nonlinear functional 0[n±{r)], given by Eq. (Fl), only up 
to the quadratic terms and neglecting (consistently under linearization) the quadratic contribution in the 
approximated averaged Euler-Lagrange equations, Eq. (F4). Although internal self-consistency (within the 
semi-grand-canonical ensemble) is achieved by using the linearized self-consistent averaged densities (|60|), 
we show in Appendix ^ that global self-consistency under linearization (between the canonical and the 
semi-grand-canonical ensembles) will also require the inclusion of the quadratic state-dependent terms of the 
self-consistent averaged densities (F5). 



G Legendre transformation at the linearized level 

In this Appendix we discuss the differences between the Legendre transformation connecting the canonical 
and the semi-grand-canonical ensembles at the level of the linearized functionals (before the functional 
minimization) and of the linearized thermodynamic potentials (after the functional minimization). It is 
shown that, in order to preserve the exactness of the Legendre transformation, quadratic contributions to 
the linearized expansion densities should be included in the former case, which are automatically included 
in the latter case. 

At the nonlinear PB level the osmotic pressuies obtained in the two distinct (canonical and semi-grand- 
canonical) ensembles are completely equivalent provided we map them using the nonlinear relation 



(n, + ns)n, = (n+(r)) (n_(r)) = (e'^W) (e 



(Gl) 



where Ug = (n_(r)) is the effective average salt concentration in the colloidal suspension. The exact (at 
the mean-field level) relation ( |G1| ) follows directly from the gauge-invariant forms of the nonlinear aver- 
age densities (|6l|). Therefore, up to quadratic order, the linearized osmotic pressures are related by the 
renormalization of the total average density of microions in the two ensembles. 



2n. 



(1 + 2,s) ne ^ Mr)) + (n_(r)) ^ !^ |l + i (l - {S,{r)) + 6 [(<53(r))]| 



(G2) 



Because of the quadratic contribution, the linearized osmotic pressures obtained in the two ensembles do not 



have the same form when they are mapped using the zero-th order renormalization (1 -|- 2s) 



In other 



words, although for the nonlinear equations the Legendre transformation between the canonical and the semi- 
grand-canonical ensembles is exact, the same does not hold for the linearized equations: one needs to use the 
approximated mapping ( |G2[ ) and expand consistently the linearized osmotic pressure up to quadratic-order 
terms. This introduces an additional source of deviations for the linearized semi-grand-canonical equations 
of state. In particular, the thermodynamically-conjugated density (in the semi-grand-canonical ensemble) 
to the chemical potential of salt particles, fig = IJ,+ + /i-, with /i± = In (nbC±)j that corresponds to the 
(effective) total average density of microions inside the colloidal suspension, n = + fi-, 
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is indeed given by the right-hand side of Eq. (G2) neglecting cubic and higher-order contributions. This 
conjugated density, however, is inconsistent (up to the quadratic order, but consistent under linearization) 
with the state-independent zero-th order Donnan densities (|6^), i.e. n ^ + fl_ = Uc/t], because of the 
presence of the quadratic contribution in Eq. ( |G2|) . 

An alternative procedure to the Legcndre transformation of the linearized Helmholtz free-energy func- 
tional ^T:)H[n±{r)] — which is done, as presented in Appendix ^, before the functional minimization with 
respect to the profiles — is to perform it after the functional minimization, at the Helmholtz free- energy 
level. Because in the canonical ensemble the expansion densities n± are known a priori, the Legendre 
transformation that maps the linearized Helmholtz free energy Fdh into the linearized semi-grand-canonical 
potential Quh can then be obtained without any further approximations for the expansion densities. For 
this purpose it is convenient to recall the definition of the total volume of the WS cell and of the nominal 
densities of counterions, polyions and salt particles. 
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It is important to stress that the nominal densities should be used, instead of the effective ones. Introducing 
the linearized Helmholtz free-energy density, /dh = Fbh/Vj where ^dh in the presence of added salt is given 
by Eq. (^), one may check that the Legendre transformation at the linearized free-energy level is indeed 
exact, since the linearized chemical potentials of salt particles and of polyions. 
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satisfy the thermodynamical identity 



^DH = ^sfJ-s + np^p ~ /] 



DH, 
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where the linearized canonical osmotic pressure entering into Eq. (G7), P^m, is given by Eq. (pO[). It 
should als o be remarked that Eq. ( |G5| ) corresponds to the truncation of the expansion of the exact nonlinear 
mapping (Gl) up to quadratic-order terms. 
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The thermodynamical relation ( G7 ) can also be viewed as defining the Legendre transformation. Instead 
of obtaining the osmotic-pressure isotherms for a constant number of salt particles inside the WS cell (canon- 
ical case, fixed s), we may consider them at fixed chemical potential of salt particles /is (semi-grand-canonical 
case), which corresponds to a system in electrochemical equilibrium with an infinite salt reservoir of bulk 
density n\,, defined by 

= + = In (n^C^C!) . (G9) 

Solution of the nonlinear equation defined by ( |G5[) and (|G9| ) yields the Legendre transformation in an implicit 
parametric form. 
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with the quadratic electrostatic-potential deviation in the canonical ensemble given by 

1 1 Znis 



Mr)) = 



(l + 2.s)^ 1 + 2s 2A|(Ki?, Ka) 

Al(Ki?, KO) 



[Ai(Ki?, K.a) - A2{kR, na)] - 4Ka (l - k^E^) - Ak^R^ 



(Gil) 



The linearized semi- grand- canonical osmotic pressure, = ^^^(0, nb), is then obtained by inserting 



the Legendre transformation s = s{(j),n\j) into the linearized canonical osmotic pressure 



pcan 
^DH 



Eq. (50). In other words, the two linearized osmotic pressures are related by P^^ifp, ^h) = Pim i'P^ ^(4'' "■b)]- 



We should remark that the linearized semi-grand-canonical osmotic pressure obtained by this procedure, 
Pq^((/>, nb), does not coincide with its counterpart ^^^(0, nb) given by Eq. @. 

The disagreement between the two distinct linearized semi-grand-canonical osmotic pressures, ^dhI*?^' ^b) 7^ 
Ppy((/), rib) — obtained by Legendre transformations performed pre- and pos-minimization of the linearized 
functional ^dh["-± (''')] with respect to the profiles — may be traced back to the neglected quadratic contri- 
bution in the self-consistent linearized averaged densities (|60|). Although the neglect of this state-dependent 
contribution in the average densities ( p^ ) is justified to ensure mathematical consistency of the truncation 
under the linearized approximation, in order to obtain global thermodynamic self-consistency (i.e., in order 
to preserve the exactness) of the Legendre transformation one needs to keep all terms of the quadratic trun- 
cation of the Euler-Lagrange equations, regardless of its inconsistency from the mathematical point of view. 



This leads to the average densities (F5), which include the quadratic state-dependent contribution {S2{r)) 
— in contrast to Eqs. (|60[), which neglect it. Compared to the linearized semi-grand-canonical potential 
flDH, Eq. (p9|), which uses the state- independent zero-th order Donnan densities (|60|) as expansion densities, 



an augmented version r^DH using the quadratic average densities (F5) will have an additional quadratic 
state-dependent contribution, 
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where the parameter t) and the (effective) Debye screening length k ^ in the suspension. 
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are now given implicitly in terms of the quadratic electrostatic-potential deviation in the semi- grand- canonical 



ensemble, which is obtained by replacing (1 -I- 2s) 



and 



k in the expression of its canonical 



counterpart, Eq. ( Gil ). These yield the globally self-consistent (i.e., that preserve the exactness of the 
Legendre transformation) linearized semi-grand-canonical equations of state, 
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where the total derivatives must take the {62{r)) parametric implicit dependence of Eqs. (G13) into account. 
The (effective) total average density of microions n, Eq. (G14), the linearized semi-grand-canonical osmotic 



pressure P^^, Eq. ( G15| ), and the linearized chemical potential of polyions /ip((/), n^,) are now fully consistent 

" " They 



with their canonical counterparts, given respectively by n = (1 + 2s)nc and Eqs. ( |50| ) and ( p6[) . 1 
are related by the Legendre transformation (1 -I- 2s) = where ij — given implicitly by Eq. ( |G13 ) 
also includes quadratic state-dependent contributions. A further Legendre transformation of the linearized 
semi-grand-canonical potential rion regains, as it should, the linearized pressure P^^, 



psgc 
-^DH 



DH 



V 



(G16) 



25 



The globally self-consistent linearized osmotic pressure P^^, Eq. ( |G15 ), which preserves the exactness 
of the Legendre transformation, leads to results qualitatively similar to those of Eq. ([to]). In particular, 
the nonmonotonic behaviour of the spinodal lines for weak screening (^ba <C 1) and the intrinsic instabil- 
ity of the low-(f> phase are still predicted by P^^, as shown in Figure ^, where we compare the spinodal 
lines associated with the two distinct semi-grand-canonical linearized osmotic-pressure definitions, given 
by Eqs. ( |70| ) and ( Gl|5| ). We should mention, however, that explicit analytical comparison in the exactly 
solvable planar casa23 does not show any improvement of the agreement between the nonlinear and lin- 
earized equations with the inclusion of the quadratic contribution to the aiterage densities. Any numerical 
indications in this direction, which were indeed observed in the planar case,c3 are purely fortuitous. In fact, 
asymptotic expansions in the weak-coupling (£s — > 0) and in the ideal-gas limit of both linearized osmotic 
pressures in the planar case, and P^^, agree with the full nonlinear PB version up to the same order. 
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Note that the nonhnear osmotic-pressure difference satisfies /3AP = n+(i?) + n-{R) — 2nb = 
4nbsinh2{i[/3^ei-^(i?)]} > 0. 

Handbook of Mathematical Functions, edited by M. Abramowitz and I. A. Stegun (Dover, New York, 
1970). 

This equivalence does not hold beyond the mean-level, because microionic correlations in the bulk salt 
reservoir will lead to different osmotic pressures in the two distinct cases. 
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Figure 1: Salt-free (s ~ 0) linearized osmotic-pressure isotherms as a function of the volume 
fraction (f) ~ {a/Kf. From top to bottom the isotherms correspond to 6* = 41, 43, 6*0111 = 
44.902477- • • (bold line), 47 and 49. In the gray region the salt-free linearized isothermal compressibility 
Xdh is negative, which would imply a thermodynamical instability that leads to a phase separation between 
two fluid phases: a low-0 (dilute) gas (G) and a high-c/i (dense) liquid (L). The black circle represents the 
salt-free critical osmotic pressure and the dashed curve defines the salt-free spinodal line in the 0x0 diagram 
(the s — Q line in Figure 
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Figure 2: Spinodal lines associated with the linearized canonical osmotic pressure Poff, in the 9 = 'iZis/a 
versus volume fraction = [a/ Kf' plane. They correspond to lines of vanishing inverse isothermal com- 
pressibility, Xrm ~ ^- gray region the linearized isothermal compressibility of the salt-free (s ~ 0) 
suspension becomes negative, leading to a coexistence iietween gas (G) and liquid (L) fluid phases. Note 
that this is in contrast to the full nonlinear treatment,E3 which always predicts positive compressibilities. 
Addition of monovalent salt reduces the unstable region by shifting the spinodal lines to higher values of 0, 
as labeled by the different curves with increasing values of s. The black circle represents the salt- free critical 
point (see main text for more details). 
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Figure 3: Spinodal lines (xdh ^ ^) associated with the hnearized semi-grand-canonical osmotic pressures 
Po'g (dashed hues), Eq. @, and P^^ (sohd hues), Eq. ( |gT5| ), in the d = 3Z£b/ a versus volume fraction 
4> — {a/Kf plane. They delimit the spurious unstable region that extends to lower values of 4>. To allow 
a comparison with the canonical case (Figure ||), we also show the salt-free critical point (black circle) and 
the salt-free (^ba = 0) unstable gray region. In the salt-free limit (^ba <C 1) the semi-grand-canonical 
spinodal line reduces to the salt-free one, although for any nonvanishing n^a eventually it will bend to 
the zero-temperature critical point at ((/icrit, ^crit) = (0,0). Contrary to the canonical case, an increase of 
the bulk reservoir density in the semi-grand-canonical case enhances the instability, as can be seen from the 
different spinodal lines with increasing KbO. A typical monotonic (non-oscillating) osmotic-pressure isotherm 
is presented in Figure 3 (dotted curve) in Ref. |^^. See also main text for additional comments. 
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